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ABSTRACT 

Context. MWC 656 has recently been established as the first observationally detected high-mass X-ray binary system containing a 
Be star and a black hole (BH). The system has been associated with a gamma-ray flaring event detected by the AGILE satellite in July 
2010 . 

Aims. Our aim is to evaluate if the MWC 656 gamma-ray emission extends to very high energy (VHE > 100 GeV) gamma rays. 

Methods. We have observed MWC 656 with the MAGIC telescopes for ~23 hours during two observation periods: between May and 
June 2012 and June 2013. During the last period, observations were performed contemporaneously with X-ray (XMM-Newton) and 
optical (STELLA) instruments. 

Results. We have not detected the MWC 656 binary system at TeV energies with the MAGIC Telescopes in either of the two 
campaigns carried out. Upper limits (ULs) to the integral flux above 300 GeV have been set, as well as differential ULs at a level of 
~ 5% of the Crab Nebula flux. The results obtained from the MAGIC observations do not support persistent emission of very high 
energy gamma rays from this system at a level of 2.4% the Crab flux. 

Key words, binaries: general - gamma rays: observations -gamma rays: binary - stars: individual (MWC 656) - X-rays: binaries - 
X-rays: individual (MWC 656) 


1. Introduction 

High-mass X-ray binaries (HMXBs) are systems composed of 
a massive star (M* > 10 M©) and a compact object, either a 
black hole (BH) or a neutron star (NS). The search for GeV 
and TeV emission from HMXBs has been the aim of exten¬ 
sive studies during the past few decades. Despite the large num¬ 
ber of observations devoted to the search, only a few of those 
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systems have been confirmed as gamma-ray emitters. A partic¬ 
ular group of five systems are reg ularly detecte d at TeV ener¬ 
gies: the gamma-ray binaries (see lDubusll2013l and references 
therein). Two other HMXBs have been the object of extensive 
searches: Cygnus X-3, which e mits in the high-energy (HE; 100 
MeV < E < 100 GeV) dom ain dEermi-LAT Collaboration et al.l 
120091 : iTavani et all 12009^ . and Cvg X-1. which has been 
repor ted to emit at H E (ISahatini et all 1201 3l iMalvshev et al.l 
1201 3l Bqdaghee^t^ 2 OI 3 I) and showed a ~ 4cr excess at 
VHE ( Albert et al. 20071) . To investigate the gamma-ray mech- 
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anisms in this type of sources, observational campaigns on 
other HMXBs have been carried out. Th e recently discovered 
object MWC 656 (iLucarelli et al.l 1201 Oh is a HMXB system 
and has been propos ed as a new gamma-ray binary candidate 
(IWilliams et al.ll2010h . _ _ 

On July 2010, AGILE (iTavani et al.l l2009ah detected a 
gamma-ray point-like source dubbed AGL J2241-1-4454 with a 
significant excess above 5 sigma, displ aying an integral flu x 
above 100 MeV of 15 x 10“^ ph cm“^ s ' (ILucarelli et al.ll201^ . 
The source was first detected during the period between 25th 
July at 01:00 UT (MJD = 55402.042) and 26th July 2010 at 
23:30 UT (MJD = 55403.979). The source is located at (/,h) = 
(100.0°,-12.2°) + 0.6°(95%stat.) + 0.1°(syst.). At the time of 
writing, no further flares from this source have been reported 
and no spectru m has been published . 

Term/-LAT (I Atwood et al.ll20()^ could not confirm the de¬ 
tection by AGILE and an analysi£] of simultaneous data from 
the same direction yielded an upper limit (UL) of 10“^ ph cm“^ 
s ' (95% confidence level, from now on CL) above 100 MeV, 
assuming a photon index T = 2. A more extended analysis 
of Termi-LAT data, including 3.5 years of data on the AGL 
J2241-1-4454 source location, also led to no evidence of HE 
gamma-ray emission. A 90% CL UL was set at t he level of 
9.4 X 10^'° ph cm“^ s ' for 3.5 years of observations (iMori et al.l 
1201 3h . 

The Be star MWC 656, also known as HD 215227, lies 
within the error bars of the AGILE best-fit source position. It 
was proposed as the optical counterpart of the ex cess claimed 
by the AGILE collaboration (IWilliarns et al.ll201^ . The system 
displays opti cal photometric modulation with a period of 60.3 7 
+ 0.04 days (IWilliams et al.ll20Tot IParedes-Fortunv et al.ll2012l) . 
Optical spectr oscopic measureme nts of MWC 656 confirmed its 
binary nature dCasares et al.ll201^ . Recent optical spectroscopic 
measurements improved the spectral classification and reduced 
the uncertainties in the spectrophotometric distance, placing the 
system at a distance of 2.6 + 0.6 kpc. These measurements also 
revealed that the compact object is a stellar-mass BH of 3.8-6.9 
solar masses, making th is the first known case of a Be/BH sys¬ 
tem (ICasares et al.ll2()14h . 

MWC 656 was also observed in radio with the Europ ean 
VLBI Network (EVN) and was not detected: [MoldonI ( |20 1 2b re¬ 
ported 3cr radio flux density ULs at 30-66 piy level. 

X-ray observations were performed by XMM-Newton 
when the source was a t an orbital phase (p = 0. 080 
dMunar-Adrover et al.l2014t) . The X-ray flux measured was com¬ 
pared with the radio ULs , resulting in a ratio compatible with 
the correlation derived in ICorbel et akl d2013l) for BH LMXBs, 
and comparable to the faintest BH LMXBs detected. A search 
for hard X-ray emission has been conducted with INTEGRAL 
dLi et al.ll2013[) with no positive detection in the 18-60 keV en¬ 
ergy band reported for a total exposure time of 2.1 Ms. 

In addition, the MAXI mission, which continouosly monitors 
the X-ray sky in the 2-20 keV band, has not detected emission 
coming from the AGL J2241-1-4454 positior0 on the same date 
as of the AGILE detection. 

In this work we present the results of the observations of 
MWC 656 carried out with the MAGIC telescopes in 2012 and 

’ http://fermisky.blogspot.com.es/2SlS/07/ 
extra-note-july-30-2SlS.html 

^ Phase 0 has been set to the maximum of optical brightness, on 
HID 2453243.3 (MJD 53242.8). With the ephemeris from ICasares et al.) 
( 120141) . the periastron passage occurs at phase 0.01 ± 0.10. 

^ http: //maxi . riken. jp/top/index. php?cid=l(§;jname= 
J2242+447#lsp 


Table 1: Observations of MWC 656 performed by MAGIC in 
2012 and 2013. 


Date 

(MJD) 

Orbital 
Phase (°) 

Zenith Angle 
Range (°) 

Time 

(hours) 

Mode 

56070 - 56078 

0.83 - 0.95 

23-50 

9.4 

mono 

56092 - 56097 

0.20 - 0.28 

22-51 

14.0 

mono 

56446 - 56448 

0.06 - 0.08 

28-45 

3.3 

stereo 


2013. X-ray and optical observations were performed during the 
2013 campaign to study the behavior of the source in a multi¬ 
wavelength context. 


2. Observations 

VHE observations of MWC 656 were carried out using the 
MAGIC telescopes, which are located at the observatory of El 
Roque de Los Muchachos (28°N, 18°W, 2200 m above the sea 
level) on the island of La Palma, Canary Islands, Spain. The 
system consists of two 17 m diameter Imaging Atmospheric 
Cherenkov Telescopes (lACTs) each one with a pixelized cam¬ 
era containing photo-multipliers, covering a field of view of 3.5°. 
The current sensitivity of the MAGIC stereoscopic system is 
0.71% + 0.02% of the Crab Nebula flux in 50 h of observation 
for energies above 250 GeV (lAleksic et alJl201^ . The spatial 
resolution at these energies is < 0.1° and the energy resolution 
is ~18%. In the case of monoscopic observations (also referred 
as mono observations) the integral sensitivity above 280 GeV 
is abo ut 1.6% of the Crab Nebula flux in 50 hours (lAliu et al.l 
l2009l) . The observations are performed using wobble mode, in 
which the telescopes point at two different symmetric regions 
situated 0.4° away from the source position to simultaneously 
evaluate the background. 

The observations of MWC 656 were performed during two 
different epochs: May-June 2012 and June 2013. The observa¬ 
tions in 2012 were performed between 23rd of May and 19th 
of June in mono mode with MAGIC-II (the MAGIC-I telescope 
was not operational) for 23.4 hours. After selecting good-quality 
data, a total of 21.3 hours remained. The 2013 observations were 
performed between June 3rd and 5th, just after the periastron 
passage (see Figure[T]i in stereo mode. The source was observed 
for a total of ~3.3h during this period. A summary of the obser¬ 
vations is shown in Table [1] 

The observation of June 4, 2013 {(p - 0.08) lasted for ~ 1 h 
and was taken almost simultaneously with a XMM-Newton ob¬ 
servation {XMM-Newton started observing right after MAGIC 
finished its observations), t he results of which are reported in 
llVlunar-Adrover et al.l (120141) . MWC 656 was also observed with 
the fiber-fed STELLA Echelle Spectrograph ( SES) of the 1.2m 
robot ic STELLA-I (ST) optical telescope (IStrassmeier et al.l 
l2004h at the Observatorio del Teide in Tenerife on the nights of 
2nd, 3rd, 5th and 8th of June 2013. The spectra cover the wave¬ 
length range 3870-8800 A with increasing inter-order gaps start¬ 
ing at 7200 A. The spectrograph provides an effective resolving 
power of R =55000. Two spectra were obtained on the nights of 
2 and 5 June and one on the nights of 3 and 8 June. The inte¬ 
gration time was set to 1800 s per spectrum while the automatic 
pipeline products were used for the extraction and calibration of 
the spectra. 
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Fig. 1; Depiction of the orbit of the MWC 656 system as seen 
from above. The optical star MWC 656 lies at the focus of the 
ellipse and the BH follows the elliptical orbit. The size of the 
star is scaled with respect to the BH’s orbit. The MAGIC, XMM- 
Newton and AGILE observations are overlaid along the orbit. 
Circles represent steps of 0.1 orbital phases while triangles mark 
the periastron and apastron phases, which are linked by a dot¬ 
ted line. Squares mark the position of the inferior and superior 
conjunctions. 

3. Data Analysis 

The MAGIC data analysis was performed using the stan- 
dard MAGIC ana lysis and reconstruction software, MARS 
dZanin et al.l l2013l) . The record ed shower images were cali - 
brated, cleaned and parametrized (lHillaslll9^lAliu et al.ll2009l) . 
The y/hadron separation (backgro und rejection) is pe rformed via 
the Random Forest (RF) method (lAlbert et al.ll2008b . The event 
direction and the energy of the primary gamma ray were recon¬ 
structed, in the case of mono observations, also by using a RF 
method. The energy of each event in the case of stereoscopic ob¬ 
servations is estim ated using look-up t ables generated by Monte- 
Carlo simulations (lAleksic et al.l201^. Upper limi t s (UL s) were 
derived using the method explained in iRolke et al.l (l2005l) with a 
CL of 95% and a systematic uncertainty of 30%, assuming dif¬ 
ferent photon indexes (F = 2.0, 2.5 and 3.0). The values obtained 
for the three spectral indexes are compatible at the 5% level. The 
results reported in this paper are for F = 2.5. 

4. Results 

No signihcant gamma-ray emission was detected from 
MWC 656 in either observational campaign. Furthermore, no 
significant signal was detected in a day-to-day analysis. 

We computed 95% conhdence level (CL) integral flux ULs 
above 300 GeV. The integral flux UL for the whole observational 
campaign of MWC 656 is F(E > 300 GeV) = 2.0 x 10 cm“^ 
s“* (-2.4% of the Crab Nebula flux at the same energy and cor¬ 
responding to a luminosity Lvhe ~ lO^^erg s“*) at 95% CL, 
assuming a power-law model with a photon index F = 2.5. 

We divided the observational periods into four different 
phase bins, using a bin width of 0.1 in phase, along with the most 
recent ephemeris (ICasares et al.ll2014l) . The phase was binned as 
follows: phases 0.8-0.9, 0.9-1.0 and 0.2-0.3 for the 2012 cam¬ 
paign, covering the orbit before the periastron passage and also 
post-periastron phases (see Figure[TJ. The 2013 campaign cov¬ 
ered the phase range 0.0-0.1, just after the periastron passage. 


Table 2: Integral flux ULs for E > 300 GeV calculated at 95% 
CL for MWC 656 for each orbital phase range. 


Mode 

Phase bin 

Integral UL 
{E >300 GeV) 

Significance 

tv/ 



(10“^^ cm“^ s“') 

(cr) 

(h) 

stereo 

0 .0-0.1 

2.0 

1.0 

3.3 

mono 

0.2-0.3 

8.7 

2.1 

4.9 

mono 

0.8-0.9 

6.5 

1.0 

11.5 

mono 

0.9-1.0 

2.5 

-1.1 

4.9 


No signal was detected and integral ULs for bins of ~0.1 in phase 
have also been computed. (See Table |2]) 

We have computed differential flux ULs from the energy 
threshold of our analysis (245 GeV) up to 6.3 TeV at 95% CL, 
with five bins per decade of energy (see Figure|2]). 

The MAGIC observations carried out on June 4, 2013 were 
performed almost simultaneously with an XMM-Newton obser¬ 
vation. The detected low X-ray flux was consistent with the 
source being in the quiescent state (defined in terms of the Ed- 
dington luminosity, when L < lO^^LEdd) during the observation 
dMunar-Adrover et al J l20 1 4h . The MAGIC integral flux UL for 
June 4 is F{E > 300GeV) < 4.9 x 10 cm“^ s“'. There is 
no specific information about the X-ray state of the binary sys¬ 
tem during the 2012 observations. Other space missions such as 
MAXI have not reported emission from MWC 656 during the 
2012-2013 campaign, which might be indicative of a quiescent 
state as well. 

Finally, the STELLA spectra, contemporaneous with the 
2013 MAGIC campaign, shows the presence of the double 
peaked He II /14686 emissi on line with an equiv alent width com¬ 
parable to that reported in ICasares et al.l (l2014li . We also detect 
other emission lines, mainly Ha, H/J and weak Eell lines wit h 
comparable strength to that measured by ICasares et alJ (l2012l) . 
Therefore, we conclude that MWC 656 is in a similar optical 
state as in past observations, the 2013 X-ray observations indi¬ 
cate a quiescent state, and hence the accretion activity should be 
very similar. 

5. Discussion and Conclusions 

We have searched for a VHE counterpart of the only known 
Be/BH binary system, MWC 656. The VHE observations per¬ 
formed by MAGIC can exclude a VHE flux based on the ex¬ 
trapolation of the emission from the AGILE detection. Assum¬ 
ing a power-law spectrum and a photon index F = 2.5, this 
emission would be ~ 4 x 10^^' TeV^' cm“^ s“' at 300 GeV 
(Lvhe ~ 2 x 10^^ erg s“^), which is well above the UL imposed 
by MAGIC. However, no flaring episodes were reported during 
the MAGIC observations, limiting the conclusions we can derive 
from the HE/VHE comparison. 

In this type of binary, several mechanisms have been 
proposed that would result in gamma-ray emission above 
tens of GeV (R emillard & McClintockI 120061 : iFendej 120061 : 
IZdziarski et al.l 20141) . Unfortunately, the lack of contempora¬ 
neous data at other wavelengths during the AGILE flare make 
conclusions on the type of emission model highly speculative. 
It is even possible that the AGILE detection was just a transient 
event of an unknown nature in the direction of the binary sys¬ 
tem but not related to it. Nevertheless, different emission levels 
could be expected depending on the state where the system is, 
i.e; quiescence or accreting state. 
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Fig. 2: SED of MWC 656 including MAGIC ULs from the 2013 
campaign together w i th sim ultaneous XMM-Newton data from 
iMunar-Adrover et ^ (|2014|) . We also include E VN radio ULs 
from Moldonl f 2012 ), the AGILE energy flux from iLucarelli et all 
(1201 ol) and the fermt-LAT UL simultaneous (green) to the AG¬ 
ILE measurement. The 3.5-year UL set by fermt-LAT to any 
persistent emission is also plotted. 


During simultaneous X-ray and VHE observations the X- 
ray luminosity of the source in the 0.3 - 5.5 keV energy range 
was Lx(0.3-5.5 keV) ^ (l-6!^:^) ^rg s^' = (3.1 + 

2.3) X 10 ®LEdd (iMunar-Adrover et al. 201^ for the estimated 
BH mass range 3.8-6.9 M© ( Casares et alJ 20141) . The low X- 
ray luminosity is characteristic of systems in quiescent states 
(dehned i n terms of the Ed dington luminosity, when Lx < 
lO^^LFHH. lPlotkin et al.ll201^ . For instance, the X-ray luminos¬ 
ity is ~5 orders of magnitude lower than the one typically ob- 
served in Cygnus X -1, which has also been observed by MAGIC 
(lAlbert et al.l2007l) . Even if we consider an increase in the X-ray 
flux consistent with a flaring state in the 2012 observations, and 
a ratio between X-rays and VHE gamma ra ys of Lx/LtrV ~ 1 0 
(similar to that observed for Cygnus X-1 in lAlbert et al.ll206^ . 
the expected VHE emission would be a factor ~ 1.5 x 10“^ of 
the Crab Nebula flux ( ~ 2 x 10^“ cm“^ s“' ) in this case, also 
well below the detectable levels for the current lACTs and even 
for the next generation of Cherenkov telescopes; the Cherenkov 
Telescope Array (CTA). The integral sensitivity of CTA is fore¬ 


3xl0-“TeV 


cm 


^ above 50 GeV for 50 h of 
Mn 1 h observation, con- 


seen to reach 

observation (~ 7 x 10^“ TeV __ 

sidering array E configuration) (lAcharva et al.l 120131) . still not 
enough to detect MWC 656 in relatively short time in this sim¬ 
ple approximation. 


Eigure |2] shows the spectral energy distribution (SED) of 
MWC 656. The absence of a detection of steady emission at 
high and very high energies implies that the X-ray emission can¬ 
not continue to increase with energy indefinitely and must turn 
over in the SED. The MAGIC differential ULs correspond only 
to June 2013 data, because for the 2012 observations X-ray in¬ 
formation is not available. We also plotted the AGILE measure¬ 
ment in the SED along with the Fenni-LAIT upper limit, obtained 
with observations performed on the same dates of the AGILE 
detection. Although the Lermf-LAT UL contradicts the AGILE 
detection, it is worth noting that the observation modes of these 
telescopes are different and that the integration time might not be 
exactly the same. Therefore, the observations are not strictly si¬ 


multaneous and Fermi-LAH might have missed short ( < 1 hour) 
gamma-ray flares from MWC 656. 
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